Spatial coherence in complex photonic and plasmonic systems 
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The concept of cross density of states characterizes the intrinsic spatial coherence of complex 
photonic or plasmonic systems, independently on the illumination conditions. Using this tool and 
the associated intrinsic coherence length, we demonstrate unambiguously the spatial squeezing of 
eigenmodes on disordered fractal metallic films, thus clarifying a basic issue in plasmonics. 
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The optical properties of nanostructured materials 
have attracted a lot of attention, due to their potential 
for light concentration and transport at subwavelength 
scales [TJ [2]. New possibilities have emerged, for the 
design of efficient sources and absorbers of visible and 
near-infrared radiation, or for optical storage and infor- 
mation processing with ultrahigh spatial density. Metal- 
lic nanostructures benefit from the excitation of surface 
plasmons that permit concentration at ultra-small length 
scales and ultra-fast time scales [3]. Periodic materials 
structured at scales on the order of the wavelength (pho- 
tonic crystals) [4 or below the wavelength (metamateri- 
als) [5] benefit from bandgap effects or artificial effective 
properties going beyond that of natural materials. Disor- 
dered media also offer the possibility to build up spatially 
localized modes (e.g. by the process of Anderson local- 
ization) [6]. Light concentration and transport at sub- 
wavelength scales encompass a broad range of processes, 
including coherent control at the nanoscale [7], enhance- 
ment of light-matter interaction in weak and strong cou- 
pling regimes [5HT2] , superradiance [13] , enhancement of 
non-radiative energy transfer [14], or light focusing be- 
yond the diffraction limit [T5UT7] . The spatial extent of 
eigenmodes is of central importance, since it characterizes 
the ability of the system to support concentrated or delo- 
calized excitations. It drives, e.g., the coherence length of 
surface plasmons [12 , 18-21], the range of non-radiative 
energy transfer [22 , 23 , or the lower limit for spatial fo- 
cusing by time reversal or phase conjugation [24H26] . The 
trade-off between localized and delocalized excitations is 
also a central issue for the understanding and the con- 
trol of the optical properties of disordered fractal metallic 
films [27 . In this Letter, we introduce the Cross Density 
Of States (CDOS) as a quantity that characterizes the 
overall spatial extent of eigenmodes, and use it to ad- 
dress the spatial localization of light on disordered frac- 
tal metallic films. We demonstrate unambiguously the 
spatial squeezing of eigenmodes close to the percolation 
threshold, thus providing a theoretical basis to clarify a 
controversial issue in plasmonics [8 , 28-30 . This also il- 
lustrates the relevance of the CDOS to characterize the 
intrinsic spatial coherence in photonics and plasmonics 
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systems. 

In order to characterize the intrinsic spatial coherence 
of complex photonic or plasmonic system at a given fre- 
quency lj, we introduce a two-point quantity p(r, r',u;) 
that we will refer to as CDOS, defined as 

p(r, r', u) = ^Im [Tr G(r, r', to)} . (1) 

In this expression, c is the speed of light in vacuum, 
G(r, r',uj) is the electric dyadic Green function that 
connects the electric field at point r to an electric- 
dipole source p at point r' through the relation E(r) = 
Pquj 2 G(r,r / ,cj)p, and Tr denotes the trace of a tensor. 

The choice of this quantity as a measure of the intrin- 
sic spatial coherence results from the observation that 
the imaginary part of the Green function at two differ- 
ent points appears in a number of situations where the 
spatial coherence of random fields (produced by random 
sources and/or a disordered medium) needs to be char- 
acterized [6j [3TM33] . The imaginary part of the Green 
function also describes the process of focusing by time 
reversal in a closed cavity [24j [25]. The precise defini- 
tion of the CDOS in Eq. ([!]) has been chosen so that it 
reduces to the Local Density Of States (LDOS) when r 
and r / coincide [2] [34] . 

The physical picture behind the CDOS is a counting 
of optical eigenmodes that connect two different points 
at a given frequency. In a network picture, the LDOS 
measures the number of channels crossing at a given 
point, whereas the CDOS measures the number of chan- 
nels connecting two points. In order to give a more rig- 
orous basis to this picture, let us consider the canoni- 
cal situation of a non-absorbing system (e.g., a nanos- 
tructured material) placed in a closed cavity. In this 
case an orthonormal discrete basis of eigenmodes can 
be defined, with eigenfrequencies uj n and eigenvectors 
e n (r). The expansion of the dyadic Green function reads 
G(r,r / ,cj) = J n c 2 e;(r')e n (r)/(^ - uj 2 ) where the su- 
perscript * stands for complex conjugate [35]. This sin- 
gular expansion can be dealt with using lim r? ^ + W/( x ~ 
xq — if])] = P[l/(x — xq)] + in S(x — xq) where P stands 
for principal value, and the CDOS defined by Eq. ([I]) can 
be rewritten as 

p(ry, W ) = ^Tr[e n (r, W )<(r>)] 8(u> - w n ) . (2) 
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This expression explicitly shows that the CDOS sums 
up all eigenmodes connecting r to r' at frequency cj, 
weighted by their strength at both points r and r'. In the 
case of an open and/or absorbing system, a frequent situ- 
ation in nanophotonics or plasmonics, the rigorous intro- 
duction of a basis of eigenmodes is more involved (an ap- 
proach has been proposed recently in the quasi-static ap- 
proximation [36]). Assuming weak leakage, quasi- modes 
can be introduced in a phenomenological way, by broad- 
ening the eigenmodes using a linewidth j n . This results 
in an expansion similar to ([2| with a Lorentzian lineshape 
replacing the delta function. The eigenmode expansion 
of the CDOS gives a rigorous basis to the network picture 
put forward above. It is important to note that, using 
Eq. 0, the correct counting of modes is automatically 
performed in the calculation of the imaginary part of the 
dyadic Green function, even when a basis of eigenmodes 
cannot be explicitly calculated or even rigorously defined. 

We shall now show that the concept of CDOS allows us 
to clarify an important issue in nanophotonics concern- 
ing light scattering and localization in disordered fractal 
metallic films. These peculiar structures exhibit optical 
properties that strongly differ from those of bulk met- 
als or ensembles of isolated nanoparticles [27 j. The in- 
terplay between surface-plasmon resonances and multi- 
ple scattering by the fractal percolation clusters leads to 
spatial concentration of light in subwavelength areas (hot 
spots) [37j [38] . The theoretical description of this phe- 
nomenon has been the subject of a controversy. Using a 
scaling theory in the quasi-static approximation, a mech- 
anism based on Anderson localization has been put for- 
ward [39 . This approach has been questioned by numer- 
ical simulations showing that at a given frequency, local- 
ized and delocalized plasmonic eigenmodes were expected 
to coexist, a situation referred to as inhomogeneous lo- 
calization [28] (localized eigenmodes have a spatial ex- 
tent smaller than the sample size, while delocalized eigen- 
modes spread over the whole sample). Measurements of 
intensity fluctuations in the near field have led to the con- 
clusion that both kind of eigenmodes should coexist [29 , 
in agreement with expectations resulting from numerical 
simulations [30], and that localized modes should domi- 
nate around the percolation threshold (but not exactly at 
percolation). More recently, measurements of near-field 
LDOS statistics have confirmed the existence of spatially 
localized modes in the regime dominated by fractal clus- 
ters (close to the percolation threshold) [8] . We shall now 
demonstrate that the use of the CDOS unambiguously 
describes this spatial squeezing of eigenmodes, providing 
a solid theoretical basis to the most recent observations. 
We introduce the intrinsic coherence length as a natural 
measure of the overall spatial extent sustained by the full 
set of eigenmodes. This clarifies a fundamental issue in 
the optics of disordered fractal films. 

The CDOS can be calculated numerically using exact 
three-dimensional simulations. We summarize the proce- 
dure that is fully described in Ref. [40 . Semi-continuous 
gold films are generated using a Kinetic Monte-Carlo 



algorithm, reproducing the geometrical features or real 
films. Typical realizations of films are shown in the top 
row in Fig. [I] (with gold in black color). To calculate 
the electric dyadic Green function G(r, r', cj), we need to 
calculate the electric field at a position r generated by a 
point electric-dipole source p located at position r' . To 
proceed, we solve numerically the Lippmann-Schwinger 
equation: 

E(r)=E (r) + ^[e( W )-l] / G (r,n, w)E(n)d 3 n 

c JV 

(3) 

where Go is the vacuum dyadic Green function, Eo(r) = 
/io^ 2 Go(r, r', u)p the incident field, e(uu) the dielectric 
function of gold and V the volume occupied by the gold 
film. The full dyadic Green function is deduced from 
E(r) = /io^ 2 G(r, r', cj)p. The computation of the LDOS 
and CDOS follows from Eq. 0. 

We show in Fig. [I] the LDOS maps (middle row) and 
CDOS maps (bottom row) computed in a plane at a 
distance z = 40 nm above two different films (shown 
in the top row) corresponding to two different regimes. 
For / = 20% (left column), the film is composed of iso- 
lated nanoparticles whereas for / = 50% (right column) 
the film is close to the percolation threshold, a regime 
in which fractal clusters dominate (multiscale resonant 
regime) [8 , 27, 40 . Before studying spatial coherence and 
the extent of eigenmodes based on the CDOS, let us sum- 
marize here the main features of the LDOS maps [23| [40] . 
For low surface fraction (left column), LDOS peaks are 
observed on top of isolated nanoparticles that are reso- 
nant at the observation wavelength. A correspondence 
between LDOS peaks and the position of one or several 
nanoparticles is easily made. For a different observa- 
tion wavelength (not shown for brevity), particles can 
switch on or off resonance and the position of the LDOS 
peaks changes, but remain attached to individual parti- 
cles. The sample behaves as a collection of individual 
nanoparticles with well identified surface plasmon res- 
onances. In the multiscale resonant regime (right col- 
umn), the LDOS structure is more complex. The con- 
nection between the film topography and the localized 
field enhancements (that are responsible for LDOS spa- 
tial fluctuations) cannot be established on a simple one- 
to-one correspondence, a known feature of fractal metal- 
lic films [H Ell [38]. 

The maps of the CDOS p(r,r f ,w) (bottom row in 
Fig. [I]) are displayed versus r for a fixed position r' (cho- 
sen at the center of the sample). Their meaning can be 
understood as follows: They display the ability of a point 
r at a given distance from the center point r' to be con- 
nected to this center point by the underlying structure 
of the optical eigenmodes. For example, a large CDOS 
(larger than the vacuum CDOS) would allow two quan- 
tum emitters at r and r' to couple efficiently. It would 
also ensure coherent (correlated) fluctuations of the light 
fields at r and r' under thermal excitation [31]. The 
CDOS also allows one to discriminate between two hot 
spots at r and r' that belong to the same eigenmode (or 
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fluctuations |5]. Note that in this experimental study, 
the inverse participation ratio was used to connect qual- 
itatively the spatial extent of eigenmodes to the variance 
of the LDOS fluctuations, permitting only a qualitative 
comparison with the curve in Fig. [2] (the inverse partic- 
ipation ratio and the intrinsic coherence length cannot 
be compared directly). The precise shape of the curve 
might also be influenced by finite-size effects inherent 
to the numerical calculations (the calculations probably 
overestimates the weight of extended eigenmodes). The 
behavior of Var(^ co h) is also instructive. Strong fluctua- 
tions are observed in the regime of isolated nanoparticles, 
a signature of the strong fluctuations of the distance be- 
tween individual particles that directly influences £ co h- 
Conversely, in the multiscale resonant regime, the reduc- 
tion of the fluctuations reinforces the assumption of a 



mechanism based on collective interactions that involve 
the sample as a whole. 

In summary, we have shown that the CDOS character- 
izes the intrinsic spatial coherence of a photonic or plas- 
monic system, independently on the illumination condi- 
tions. Using this concept, we have demonstrated unam- 
biguously the spatial squeezing of plasmonic eigenmodes 
on disordered fractal metallic films close to the percola- 
tion threshold. This clarifies a basic issue in plasmonics 
concerning the description of the optical properties of 
these films. This illustrates the relevance of the CDOS 
in the study of spatial coherence in photonics and plas- 
monics systems, and more generally in wave physics. 

We acknowledge Y. De Wilde, M. Kociak, V. Krach- 
malnicoff and R. Sapienza for stimulating discussions. 
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